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Classical models for dark energy can exhibit a variety of singularities, many of which occur for 
scale factors much bigger than the Planck length. We address here the issue whether some of these 
singularities, the big freeze and the big demurrage, can be avoided in quantum cosmology. We use 
the framework of quantum geometrodynamics. We restrict our attention to a class of models whose 

■ matter content can be described by a generalized Chaplygin gas and be represented by a scalar field 
' with an appropriate potential. Employing the DeWitt criterium that the wave function be zero at 
, the classical singularity, we show that a class of solutions to the Wheeler-DeWitt equation fulfilling 

■ this condition can be found. These solutions thus avoid the classical singularity. We discuss the 
. ' reasons for the remaining ambiguity in fixing the solution. 

1—5 

in 

5^ ' Understanding the observed acceleration of the universe rnay affect our understanding of gravity as well as enlarge 
bJO, the framework of particle physics, see, for example, [U, 0, S0] for some reviews on the state of art. The presence of 
precision observations is particularly promising in this respect 
' Even though a cosmological constant is the simplest way to explain phenomenologically the late-time acceleration 
^ I of our Universe, it is not entirely appealing from a theoretical point of view. The reason is the mismatch between the 
required observational value and the expected one from theoretical grounds . This has induced a whole manufacture 
' of theoretical model building, aiming at an explanation of the recent speed up of the Universe. Most of these models 
\l invoke a dark energy component with negative pressure or a modification of gravity on large scale Q which, by 
, weakening the gravitational interaction, achieve the accelerated expansion. More exotic explanations invoke ideas of 
• a multiverse and may therefore have to use the anthropic principle or isotropic and inhomogeneous cosmologies which 
' violate the cosmological principle [§, [13, [ll| . 

In practical applications and for a homogeneous and isotropic universe, which is a good approximation for our 
Universe on large scales, whatever the entity responsible for the recent acceleration may be, it can be described 
effectively through an equation of state parameter, w. This parameter w is just the ratio between the pressure and 
the energy density of the unknown entity (sometimes designated as "dark energy" ) and may, of course, be timc- 
^ j , dependent. The value of w is observationally very close to w = —1, that is, the equation of state for a cosmological 
' constant. However, if w is larger or smaller than w = —1, the future of the Universe will be dramatically different. 
Let us mention that indications of a decaying dark energy have recently been reported [l^ . 

Our main interest here lies not so much in the observational significance of such models but in their relevance for 
understanding the quantization of gravity. The search for a consistent theory of quantum gravity is among the main 
open problems in theoretical physics [l^ . One aspect is the fate of the singularities which are prevalent in the classical 
theory of general relativity. The hope is that a consistent quantum theory of gravity is free of such singularities. This 
aspect can most easily be investigated in the framework of quantum cosmology: the application of quantum theory to 
the Universe as a whole [H, Because of our restriction to homogeneous models, the corresponding framework of 
quantum cosmology does not encounter the mathematical problems of the full theory. We also emphasize that quantum 
cosmology is a non-perturbative framework and thus different from approaches employing "quantum corrections" to 
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the classical theory-'^. Within this context, some of the dark-energy models are particularly suitable. Furthermore, if 
one of the models applies to our real Universe, this would provide us with important insights into its past and future, 
which could be a singularity-free and timeless quantum world. 

Indeed, during the last years, it has been shown that some dark-energy scenarios with w < —1, dubbed phantom 
energy models [l7[, can induce a new type of singularity called big-rip singularity. This is a singularity where both 
density and pressure diverge and which is attained for a universe that expands to infinity in a finite time. The 
quantum version of such a situation was extensively investigated in [l^ (see also [l9l . [20j ) . It was found there that the 
semiclassical approximation necessarily breaks down for large-enough scale factors, and so the singularity theorems 
no longer apply. 

It was soon realized that the big-rip sin gula rity is not the unique singularity related to dark energy and that different 



types of singularities could show up [2ll . |22| . |23| | . One of them is the big- brake singularity, which is obtained for a 
universe in the future and which is characterized by an infinite deceleration; the universe comes to an abrupt halt. 
The quantum cosmology of such a model was discussed at length in [23|. It was found there that, given reasonable 
assumptions, the wave function vanishes in the region of the classical singularity, which we can safely interpret as 
singularity avoidance. From theoretical investigations such classical singularities are well known [25|, but in the 
context here they occur within models possessing observational relevance. 

In this work we are interested in the generalization of these quantum cosmological results to a broader class of 
models. We shall find that the singularity avoidance is there effective, too. We are mainly concerned with two 
types of singularities known as the big- freeze^ and h\g- demarrage singularity, respectively [231. Let us give a brief 
characterization of them: 

• The big- freeze singularity (or typc-III singularity in the nomenclature of [13]) takes place at finite scale factor 
and finite cosmic time in a (flat) Friedmann-Lemaitrc-Robcrtson- Walker (FLRW) universe. At this singularity, 
both the Hubble rate and its cosmic derivative blow up. (This is not the case for a big-brake singularity.) Such 
a singularity can be induced by a generalized Chaplygin gas (GCG) [2§|. The GCG is a perfect fluid which 
satisfies the following polytropic equation of state [30| : 



where A and (3 are constants. This equation of state was introduced in cosmology with the intention to unify 
the dark sectors of the Universe, that is, dark matter and dark energy [s^l. How this is achieved can be seen 
from the conservation of the energy-momentum tensor of such a fluid in a homogeneous and isotropic universe. 
It implies 



/ B \^ 

P=\^+7^)) ' (2) 

where i? is a constant. Therefore, if A, B, and 1 + (3 are positive, the energy density p interpolates between dust 
energy density for small scale factor and a constant energy density for large scale factor. However, it is possible 
that the GCG can exclusively correspond to dark energy in a FLRW universe. 

The behavior of a GCG can be quite different if the parameters A, B and \ + (3 are not all positive. In particular, 
it can induce different sorts of singularities [3l|. Moreover, it was shown that a phantom GCG, that is, a fluid 
which satisfies the polytropic equation of state ([T]) for p > and p + p < 0, can induce a future big- freeze 
singularity [1^. It was also realized that even a GCG fulfilling the null, strong, and weak energy conditions can 
lead to a big-freeze singularity [3l| in the past. 

• In Ref. it was also shown that a sudden singularity (or type-II singularity in the nomenclature of [l^l) can 
be induced by a GCG. A sudden singularity is characterized by the fact that the Hubble rate is finite, while its 



^ For reviews on cosmological singularities and string theory, cf. [Tsl . [Tsll and references therein. 

^ It has been recently shown that a big-freeze singularity might be simply an indication that a brane is about to change from Lorcntzian 
to Euclidean signature l26l . even though the brane and the bulk remain fully regular everywhere. This kind of behavior can happen 
in some braneworld models where the bulk is always Lorentzian and does not change its signature. If the brane and the bulk change 
simultaneously their signature like in the models discussed in [27L |28| | , observers at the brane would not perceive to go into such a 
big-freeze singularity. 
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cosmic derivative blows up at finite scale factor. If the GCG fulfils the null, strong, and weak energy conditions, 
then the sudden singularity corresponds to a big-brake singularity (23 | which takes place in the future when 
the universal expansion is stopped by an infinite deceleration. On the other hand, if the GCG corresponds to a 
phantom fluid, the sudden singularity takes place in the past. In analogy with the terminology of the big-brake 
singularity we call this kind of sudden singularity a big- demarrage singularity because the universe starts its 
expansion with an infinite acceleration. Since the big-brake singularity has been studied in (23 | , we shall restrict 
our attention here to the phantom model exhibiting a big- demarrage singularity. 

An interesting feature of these singularities is the fact that they can occur in a macroscopic universe, that is, for 
large values of the scale factor. A pertinent question is then whether quantum gravitational effects can resolve these 
singularities. Would that mean that there could be quantum effects in the macroscopic universe? How could we 
expect these singularities to be resolved through quantum gravity? This is the major motivation standing behind 
the quantum analysis of the big- freeze and big- demarrage singularities. We shall carry out our quantization in the 
geometrodynamical framework, using the three-metric and its conjugate momentum as fundamental variables. The 
governing equation in this framework is the Whceler-DcWitt equation. One can invoke various reasons why this 
framework is appropriate for studying the question at hand'^, [37j . Perhaps the most compelling one is the fact that 
the Wheeler-DeWitt equation is the wave equation which straightforwardly leads to the Einstein equations in the 
semiclassical limit. 

Our paper is organized as follows. In Sec. [iTl we review the classical cosmology with a generalized Chaplygin 
gas leading to a past/future big- freeze singularity or a past sudden singularity. To be able to study the quantum 
behavior, the GCG has to be mimicked by a fundamental field because also the matter part should have its own 
degrees of freedom. The representation of the GCG in terms of minimally coupled scalar fields will be given in 
Sec. mil The remainder of this article is concerned with the quantum properties of these scenarios. In Sec. IIVI we 
obtain solutions for the quantum cosmological model with phantom and non-phantom GCG evolving to and from a 
big-freeze singularity, respectively. We investigate whether this singularity can be quantum gravitationally avoided. 
The big- demarrage singularity is also addressed at the quantum cosmological level. The (general) results are then 
discussed for different types of boundary conditions in Sec. |Vl We present explicit results for the quantum states that 
avoid the singularities and discuss further consequences. Sec. |Vl] gives our conclusions and outlook. In the appendix, 
we present a justification of the gravitational wave function approximation used in the paper. 



II. THE CLASSICAL BIG-FREEZE AND BIG DEMARRAGE SINGULARITIES INDUCED BY A 

GENERALIZED CHAPLYGIN GAS 



A. The big-freeze singularity without phantom matter 



We review here a particular case of the plain generalized Chaplygin gas [3l|, that is, a fiuid which satisfies the 
polytropic equation of state ([T]) and fulfils the null, strong, and weak energy conditions. Such an equation of state 
may induce a big-freeze singularity in the past. This is the case^ when A < 0, B > and 1 + /3 < 0. Then the energy 
density can be written as 

1 

(3) 
where 



-1 



3(l + /3) 



3(1 + 0) 

(4) 

denotes the minimal scale factor, which is the value where the singularity occurs. We consider a spatially flat 
homogeneous and isotropic universe filled with this sort of fluid. Then, at the minimum scale factor amin, the energy 



^ For a discussion on singularity avoidance within loop quantum cosmology, cf. [32l |: for an earlier discussion where geometrodynamical 

elements were present, cf. l33l |34|. Issl. |36|1 . 
* The conclusion reached in this section also holds for yl>0, -B<0, l + /3<0 and 1-1-/3 = l/{2n), with n some negative integer number. 

For simplicity, we shall disregard this case. 
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density blows up and so docs the Hubble rate. Similarly the pressure, which reads 

,3(1+/3)1 



p = \A\ 



-1 



(5) 



diverges when the scale factor approaches its minimum value a^in- Notice that the pressure is positive and therefore 
a FLRW universe filled with this fluid would never accelerate. This particular choice of parameters can thus not 
describe the acceleration of the current Universe. In fact, the deceleration parameter 



(6) 



is always positive. Figure [T] displays w = p/p. On the other hand, the Raychaudhuri equation implies that at a^ia 
the cosmic-time derivative of the Hubble rate also diverges. In order to show that the event that takes place at amin 
corresponds to a past big-freeze singularity, it remains to be proven that from a given finite scale factor the cosmic 
time elapsed since the singularity took place is finite. This can be done by integrating the Friedmann equation. The 
cosmic time in terms of the scale factor reads Isill 



kV3 1 + 2/3 



3(1+13) 



- 1 



1 + 2/3 
2(1 + 13) 



3(l+/3) 



F\ 1,1; 



2(1 + 



;i 



3(l+/3)^ 



(7) 



where = 8ttG, and F(5, c; d; e) is a hypergeometric function, see. for example, [ssj . In the previous expression, tmin 
corresponds to the cosmic time when a universe filled by this sort of generalized Chaplygin gas would emerge from a 



past big-freeze singularity at a 



We thus find here a singularity at finite scale factor, amin, and finite cosmic 



time, train, where both, energy density and pressure, blow up, as the Hubble parameter and its cosmic derivative do 
as well. It can easily be checked that the cosmic time elapsed since the universe emerged from the past big-freeze 
singularity until it has a given size (at a given t) is finite^, that is, t — tmin is bounded for any finite scale factor. Before 
concluding this section, we notice that although the generalized Chaplygin gas analyzed here fulfils the null, strong 
and weak energy conditions, it violates the dominant energy condition for scale factors smaller than Odom, where 



= 2 a 



(8) 



At 

fldom the pressure equals the energy density and for smaller scale factors p < p. This situation is schematically 
shown in Figure [TJ 




FIG. 1: Plot of the ratio of pressure and energy density, w, for the generalized Chaplygin gas analyzed in this section as a 
function of the scale factor [Sj]. The dominant energy condition is not fulfilled for a scale factor smaller than adom- 



^ A hypergeometric series F(fe, c; d; e), also called a hypergeometric function, converges at any value e such that |e| < 1 whenever b+c—d < 0. 
However, ifO<b + c — c!<l the series does not converge at e = 1. In addition, if 1 < b + c ~ d, the hypergeometric function blows up 
at |e| = 1 m. 
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B. The big-freeze singularity with a phantom GCG 



The big- freeze singularity discussed above may also take place in the future at a scale factor a = Omax- This is 
possible if the GCG exhibits a phantom behavior, that is, if it satisfies p + p < [1^. Eventually, the pressure will be 
so negative that the universe does not only accelerate but is even super-accelerating, that is, the Hubble rate grows as 
the universe expands. This is certainly a model capable of describing dark energy. When the scale factor approaches 
flmaxj the energy density as well as the pressure diverge in a finite future cosmic time |3l| . This phantom GCG induces 
a big- freeze singularity, although in this case the event happens in the future instead of, as above, in the past p9l.[3ll|. 
Such a future big-freeze occurs for A > 0, i? < and 1 + j3 <{). Being more precise, the energy density in this case 
reads 



while the pressure is 



where 



p = -At^ 



n \3(l+/3) 



I ^max \ 
\ a J 



a s3(l+/3)' 



I ^max A 
V a / 



3(1 + 3) 



(9) 



(10) 



(11) 



is the maximal value of the scale factor. 



C. The big-demarrage singularity with a phantom GCG 



The GCG is also known to induce future or past sudden singularities |24 
distinguished: 



l3ll |40|. The following cases can be 



1. If the GCG fulfils the null, strong, and weak energy conditions with A < 0, B > and /3 > 0, a universe 
filled with this sort of gas would face a future sudden singularity, which is just the big-brake singularity already 
mentioned above. For (3=1 this fluid was named an anti-Chaplygin gas in [23| because the case with /3 = 1 
and A > is usually called Chaplygin case. Since this model was discussed at length in [2^, we shall disregard 
it here. 

2. If the GCG is a phantom fluid with A> B <{) and /3 > 0, a universe flUed with this fluid would face a past 
sudden singularity, that is, the classical evolution starts at a finite scale factor and energy density while the 
pressure of the fluid diverges. Wc name this past sudden singularity a big- demarrage singularity as the universe 
would start its classical evolution with an extremely large acceleration due to the very negative pressure of the 
fluid. This phantom GCG corresponds to an example where a future big-rip singularity can be avoided even for 
phantom matter (soj as in this case the universe would asymptotically approach de Sitter space in the future. 
Here we are rather interested in the early behavior of a universe close to the big- demarrage singularity. 

Let us add a few more comments on the second case stated above. The energy density and pressure of the phantom 
GCG can be expressed as in ([9]) and (|10p after substituting Omax by a^nin where the initial scale factor amin can be 
expressed as in (j?]). The phantom nature of this fluid induces a super-acceleration of the universe which in particular 
implies that a universe filled with this fluid would be accelerating. 



III. THE CLASSICAL BIG-FREEZE AND BIG-DEMARRAGE SINGULARITIES INDUCED BY 

SCALAR FIELDS 



A perfect fluid is an effective description of matter. As such it is valid only on certain (large) scales. At the 
quantum level, we may need a more fundamental description of matter. We herein choose the simplest possible one: a 
minimally coupled scalar field. In the following subsections, we formulate the standard and phantom GCG in terms of 
standard and phantom scalar fields, respectively. This step is an important preparation for the quantum part because 
the wave function is defined on configurations space, that is, it will depend on the scale factor and the scalar field. 



FIG. 2: Plot of the scalar field versus the logarithmic scale factor a — ln(a/ao), where ao corresponds to the location of the 
singularity in the respective models. This plot corresponds to a past big freeze singularity where ao — amin (see Eq. (|14|l '). 



A. The big-freeze singularity driven by a standard canonical scalar field 



We first show how the GCG given in Subsection III Al can be described by a standard minimally coupled scalar field 
whose energy density and pressure in a homogeneous and isotropic universe read 



(12) 



The dot corresponds to the derivative with respect to cosmic time. Then by imposing that and satisfy the 
equation of state H]), the kinetical energy density and the scalar field potential evolve with the scale factor as 



— , Via) = -\A\rT, 

1+/3 



'a„i„\3(l+/3) 



_J_ 

1+/3 



(13) 



Therefore, the scalar field changes with the scale factor as (cf. Figure [2]) 



/On 
V ( 



(1+0) 



(14) 



where 0min corresponds to the value of the scalar field at a„iin, where the singularity occurs; we shall set 4>imn ~ 
simplicity. Finally, by combining (fT3|) and ([T4l) . the scalar field potential reads 



sinhTh (^^|l + /3||0| ] ^ 

V 2 / sinhT+7? (y^Kll + (3\\(j)\ 



(15) 



where Vi — \A\ 1+*^ /2. Notice that at the minimum scale factor or at = the potential is negative and divergent. 
In fact, the potential can be approximated in this region by^ 



2/3 



(16) 



At the lowest order, the first term in HISI I does not contribute to the approximation made in I I16II because —2/3/(1 + /3) < 2/(1 + /3) < 0. 
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FIG. 3: Plot of the potential (|15p as a function of the scalar field for the value (3 = —3. The star denotes the location of the 
singularity, KE denotes the kinetic energy of the scalar field and the arrows denote the trajectory of the scalar field. 



We recall that 1 + /? < and therefore —2(3/ (1 + /3) < —2; that is, V{4)) is at the big-freeze singularity more singular 
than an inverse square potential. This is crucial for the discussion of the quantum theory below. On the other hand, 
for large values of the scale factor or at large values of the scalar field, we have 

~ 2"TT7Jl4exp(-V3K|0|) . (17) 



The general behavior of the potential is shown in Figure [S] The scalar field starts with an infinite kinetic energy 
(at the past big-freeze singularity) climbing up the potential until it reaches its top and then starts rolling down the 
potential hill. When the potential vanishes, the scale factor is equal to adom- Hence, the dominant energy density is 
violated when the scalar field takes a value such that the potential becomes negative. 



B. The big-freeze singularity driven by a phantom scalar field 

Similarly to the case just considered, the GCG discussed in Subsection IIIBI can also be described by a minimally 
coupled scalar field, although in this case the phantom nature of the GCG implies that the scalar field does not have 
the standard kinetic term; therefore, its energy density and pressure read 

P^ = -\^'' + V{cj,) , p^ = -l4,^-V{4>) (18) 
for a FLRW universe. Equating the previous quantities to ^ and pU]) . we obtain 



1 - 



\3(l+/3) 



\3(l+/3)" 



V{a) ^ -Aw?. 



1 - 



\3(l-H/3) 



^3(l+/3)" 



(19) 



Consequently, 



= ± 



K^/3 1 + /? 



arccos 



\ a ) 



3(l + <i) 



(20) 



In Figure |4] we show how the scalar field is correlated with the scale factor. From the last equation one can read 
off that the scalar field vanishes when the scale factor reaches its maximal classically allowed value. By using the 
relations p9p and (P(7)l we find the following expression for the potential: 



1 



sinT+7? |^%^k|1-H/3||0| 



■ smT+73 



■^/3 



k|1 



(21) 
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FIG. 4: Plot of the scalar fields versus the logarithmic scale factor a = ln(a/ao), where ag corresponds to the location of the 
singularity in the respective models. The plot corresponds to a future big freeze/(past)big demarrage singularity (see Eq. (|20|l 'l. 
In the scenario with a future big freeze or a (past) big demarrage singularity, the quantity —3/2(1 + (3)a is always negative. 




KE = 



04 0^6 0^8 i 12 1^4 1^6 



FIG. 5: Plot of the potential defined in (|2ip as a function of the scalar field for the value /3 = —3. The star denotes the location 
of the singularity, KE denotes the kinetic energy of the scalar field and the arrows indicate the trajectory of the scalar field. 



where F_i = /2 and < (%/3/2)k|1 + < 7r/2 (see FigureE]). 

The big-freeze singularity at a = flmax is now located at = where the scalar field potential can be approximated 




2/3 

■ 1 + /3 



V{^)c^V^^\^K\\ + m\ ■ (22) 



This potential (also the exact potential (HH)) can be deduced from the expression (fTB|) (resp. P^ ) by rotating cj) icj) 
and taking into account that A changes sign. Notice that this implies that there is a change of sign of the potential 
(|16p after performing such an analytical continuation. 

In general, the potential ([22|) corresponds to a singular potential. However, for large enough /? (and again for 
(/) ^ 0), V{(l)) behaves effectively as an inverse square potential; that is. 



Vi^)^V^A^K\l + Pm\ , (23) 




^ Again, the second term in I I21I I does not contribute to the potential II22I I because —2/3/(1 + /3) < 2/(1 + /3) < 0. 
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FIG. 6: Plot of the potential defined in (|2ip as a function of the scalar field for the value /? = 3. The star denotes the location 
of the singularity, KE denotes the kinetic energy of the scalar field and the arrows indicate the trajectory of the scalar field. 

The scalar field starts its cosmological evolution with a vanishing kinetic energy density at -\/3/2k|1 + (}\\4>\ = 7r/2, 
then it climbs up through the potential reaching an infinite kinetic energy when the classical universe reaches its 
maximum size, that is, when ^ [see Figure [5]. This is a usual feature of phantom scalar fields. Notice as well 
that the standard canonical scalar field with the potential also climbs through its potential but in this case the 
reason behind this strange behavior is the initial infinite kinetic energy of the scalar field when (p ^ Q. 

C. The hig- demarrage singularity driven by a phantom scalar field 

Similar to the above cases one can also show that the phantom GCG leading to a big- demarrage singularity can 
be mimicked by a minimally coupled phantom scalar field (f), where 4> and the scalar field potential V{(t>) arc given, 
respectively, by Eqs ([20)1 and ([21]) after substituting Omax by Qmin, see also Figures 2] and [51 The expressions given in 
(|19p are also valid after exchanging amax by amin. The scalar field starts moving away from the singularity (located 
at = or a = amin) by rolling down the potential with an infinite kinetic energy. The phantom nature of the scalar 
field implies that (f> loses its kinetic energy as it rolls down the potential, cf. Figure [6| In fact, when the scalar field 
reaches the tail of the potential, that is, V^/2k\1 + /3||(/)| 7r/2, its kinetic energy is approaching zero, which is in 
agreement with the asymptotic de Sitter-type behavior of a universe filled by this phantom GCG. 

Close to the singularity {(p 0), V{(p) can be approximated by ([2^ . Because /? > 0, this potential does not 
correspond to a singular potential, although for 1 ^ /? it behaves as an inverse squarcpotential. Notice as well that 
for /3 = 1 the potential behaves as a Coulomb potential like its non-phantom version [2J] . 

IV. WHEELER-DEWITT EQUATION AND QUANTUM STATES 

In this section, we shall describe the quantization of the classical scenarios discussed above. This will be carried 
out in the quantum geometrodynamical framework. Its central equation is the Wheeler-DeWitt equation, depending 
on the configuration space variables (a, 4>). From the solutions obtained, we shall retrieve information concerning the 
above GCG models with regard to their quantum behavior at the classical singularity*. A discussion of the influence 
of different boundary conditions on the wave function is presented in the next section. 

Before proceeding to more technical aspects, let us mention that even though singularity avoidance is a major 
touchstone for any quantum gravity theory, no consensus regarding the criteria that may account for such an avoidance 
exists. This is mainly due to the fact that the criteria that one admits depend strongly on how one interprets the wave 
function of the universe. Divergent opinions on that interpretation arise because quantum gravity is itself a timeless 



The generalized Chaplygin has been previously analyzed from a quantum point of view in different setups (see Ref. |41| '). 
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theory [l3l |. In this respect it differs from any other quantum field theory, and there is not yet any agreement about 
the appropriate interpretational framework. 

One possible line to establish a reasonable research criterium is to regard the quantum-cosmological wave function 
as the fundamental entity from which our spacetime can be derived in an appropriate limit. The recovery of spacetime 
can be expected to occur only in special regions of configuration space. The derivation of the semiclassical limit is 
performed by a Born-Oppenheimer type of approximation scheme with decoherence as an essential ingredient, cf. 
[13| and the references therein. The gravitational part of the wave function is usually taken to be of a WKB form, 
which means that narrow wave packets around classical trajectories would not spread. But wave packets which are 
initially peaked around classical trajectories may not remain so along the entire trajectory. Such a dispersion is 
unavoidable in a quantum universe whose classical version exhibits recoUapse (4^ . Thus, a spreading of the wave 
packet signals a breakdown of the semiclassical approximation; one can then no longer associate with the wave function 
a classical spacetime as an approximate concept. Such a spreading necessarily occurs when approaching the reg ion 
of the classical big-rip singularity in quantum phantom cosmology, a phenomenon which we interpreted in [l8| as 
singularity avoidance. 

Another sufficient (but by no means necessary) criterium for singularity avoidance is the vanishing of the wave 
function at the classical singularity; there is then no possibility for such a singularity to occur in any limit. Vanishing 
of the wave function as a criterium for singularity avoidance was first suggested by Bryce DeWitt in his pioneering 
paper (43| . Singularity avoidance in this sense occurred for the big-brake singularity [23 |. Vanishing of the wave 
function will also be in the present section the appropriate criterium for singularity avoidance. 

Let us, then, turn to a detailed analysis of the quantum versions for the above discussed classical models. The wave 
function satisfies the Wheeler-DeWitt equation, which with the Laplace-Beltrami factor ordering reads 



ft2 
2 



6 da^ 



^' (a, 



-f a;je6"y(,/,)vl/(a,0) =0, 



(24) 



where F(0) is given in for the standard scalar field model and in (pij) for both phantom scalar field models. We 

have introduced the new variable a := In and assume that ag corresponds to the location of the singularity in 

the respective models. In the following, we shall use d := ^ instead of a such that ag = 1 holds. For simplicity, we 

shall drop the tilde. Recall that = 8itG. We have introduced the parameter i in order to distinguish between the 
phantom and non-phantom scalar field. For the phantom scalar field we have £ — whereas we have £ — I for 
ordinary scalar-field matter. Note that a fundamental length scale is necessary in the Wheeler-DeWitt equation to 
yield correct dimensions. To solve this equation, we make the ansatz 



where A: is a priori not restricted to real values. Furthermore, we require ipk to satisfy 



-I 



2 902 



+ ayVi(f>)cpk 



Ek{a)ipk 



(25) 



(26) 



Such a Born-Oppenheimer- type of ansatz was first used in quantum cosmology in [421 the study of wave packets. 
It assumes the approximate validity of a quasi-separability between a and (p, that is, the matter part ipk depends only 
adiabatically on the scale factor. We discuss the validity of the Born-Oppenheimer approximation in Appendix A. 

As the potentials and (|2ip are rather complicated, we solve the Wheeler-DeWitt equation for certain ranges 
of 101 and approximate the respective potential there. For the study of singularity avoidance the region of primary 
interest is \4)\ ^ 1. This corresponds in each model to the vicinity of the singularity. In this region, the potential 
is approximated by p6|) for the ordinary scalar field and by ((22)) in the case of the phantom field. Introducing the 
notation 



we find that (pk has to satisfy 



:= a^e^" 



^14101" 
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Ek{a)ipk 



(27) 



(28) 
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Defining fc^ 



and Va 



, wc finally arrive at 



213 



(29) 



where ' denotes a derivative^ with respect to (/>. We recognize that this equation is formally the same as the radial 

part of the stationary Schrodinger equation for an attractive potential of inverse power V ^ r where |0| plays 

the role of the radial coordinate r, and the angular momentum vanishes. 

Equation (|29|1 is the central equation for the following discussion. Because it formally resembles the Schrodinger 

equation, we are able to make use of results encountered in quantum mechanics. Potentials of the type V ^ r 
are there called singular [4^ . IZsj . Any potential that approaches (plus or minus) infinity faster than r~^ for r — > 
belongs to this class. For an attractive r~^-potential there exists a transitional case: if the coupling is more negative 
than a critical value, the potential is singular, otherwise regular. 

Analytical solutions for polynomial singular potentials are known for the inverse square, inverse fourth-power, and 
inverse sixth-power potentials. The inverse square potential is realized for \(3\ 3> 1, where (3 is chosen such that 
|1 -|- (3\\(t>\ is still small. The inverse fourth-power potential corresponds to /3 = —2, whereas the inverse sixth-power 
potential is realized for /3 = — |. 

In our paper we shall focus on the case \f3\ ^ 1 for the standard scalar field as well as the two phantom field models. 
We thus deal with the case of the inverse-square potential with 
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2age 



> . 



The r~^-potential in quantum mechanics was discussed, for example, in [3, [4^, [4^, cf. also [471] and the references 
therein. It became recently of interest in studying the polymer quantization which is motivated by loop quantum 
gravity [47j . 

Note that for the models with big- freeze singularity we then have /3 ^ — 1, whereas for the model with big-dcmarrage 
singularity we have /3 ^ 1. We expect that this case is sufficiently generic to accommodate also the features of other 
singular potentials. As is known from quantum mechanics (4^. [4q. the most important issue for singular potentials 
is the fact that square integrability does no longer suffice to select a unique class of states and that the spectrum 
therefore contains an ambiguity. This ambiguity points to new physics at short scales and must be fixed by either 
experiment or knowledge from a more fundamental theory. It is the fortunate property of the Coulomb potential that 
its ensuing states do not contain such an ambiguity (see, however, [ill). Thus, for the least singular potential, which 
is realized for \f3\ ^ 1, we have to solve the equation 



Ik' 



(30) 



Note that phantom and scalar matter have to obey the same quantum equation where the realm of positive energy 
for the ordinary scalar field fc^ > corresponds to the realm of negative energy for the case of the phantom field, 
fc^ < 0, cf. Eq. (|29p. The general solution to this equation is given by [49| 



<Pfc(a, W = VW\\ ci3,iV£k\^\) + C2Y,(Vik\^\) 



(31) 



J — Va, so the index is a function of a. There are four cases to distinguish: k can be real or imaginary. 



where v := 

depending on whether the energy entering fc^ is positive or negative. Furthermore, v can be real or imaginary 
depending on the parameters /3, A, and the value of a. Note that ipk satisfies 

rb 







^.(&)^(5) 



(32) 



In fact, there is one such equation for positive and one for negative (j>. Both equations are (due to the modulus-dependence of the 
potential) identical. Matching is carried out through the conditions \I'-|-(0) = ^^-(O) and = — \I''_(0) where "l/^ refer to the 

positive/negative (/(-solutions, respectively. These conditions imply that the constants in both solutions coincide. Wc will therefore just 
refer to the modulus-dependent equation I I29I I and its solution, keeping in mind that we have the same solution for positive and negative 

. 1 1 

Note that Vi = - — '-^ and V-\ = — ^ — i ^° both are positive and coincide numerically. 
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where 5 = oo for the standard scalar field and h = (hi, := ^—r-r — ^ for the phantom scalar field. (This equation is 

^ V3k|1+/3| 

analogous to Equation (11) in [43].) We have fixed the range of definition for <j) in the phantom case from to 0* 
because all classical solutions are restricted to this range. That this also holds in the quantum theory is, of course, 
an assumption, similar to treating, say, a particle on the half-line only. 

The lower bound in the integral ([32]) does not contribute as (pk vanishes at the origin (see below). Thus the matter 
dependent part of the wave function is not necessarily orthogonal. But it is if the right-hand side of ([5^ vanishes. 
This is the case if 6 = cx) and v is not an integer, as happens for the standard scalar field. The reason is the following: 
(i) The first Bessel function 3^, satisfies an orthogonality relation and (ii) the 3^{Vik\(f)\) and 3-.^{VIk\(f)\) are linearly 
independent solutions of (|30p (except when i/ is an integer). Therefore, as the solution (pk can then be expressed 
exclusively in terms of the first kind of Bessel functions, ipk is orthogonal if either of the two coefficients in front of 
3i/{Vik\(l)\) or J_i/(-\/^fc|0|) is zero. Otherwise, that is, if v is an integer, additional conditions have to be imposed to 
get orthogonal eigenfunctions, as it is the case for the phantom scalar field. 

Different types of solutions for the three models discussed in our paper will differ by the boundary conditions that 
we may choose to impose (see next section). This choice is determined by the classical trajectory. For ordinary 
matter, the classical trajectory has a minimum (see Figure 2), whereas for phantom matter with /? -I- 1 < it reaches 
its maximum at the classical singularity (see Figure 3). For the second phantom model with /3 > 0, the singularity 
lies again at a minimum of the classical trajectory. We will see how this difference in the classical model influences 
the quantum behavior through the boundary condition. Moreover, the phantom field is restricted to a finite range. 

Independent of boundary conditions, the singularity occurs in all three models at </> = and a = 0. For a = 0, 



ly ^ 1^0 y J — Va=o- There are three cases to distinguish, i — Va=o > 0, ^ — V^q^o < and the transitional case 
J — Va=o = . In the first and third case, i^q is real, whereas in the second case it is purely imaginary. 

• i - Fa=0 > 0. 

Then < ly < ^. In the quantum mechanical analogy, this corresponds to the case of a regular potential. We 
get from the behavior of the first Bessel function for small argument, 

J.W-(I)'^^^, .^-1,-2,-3..., (33) 

cf. [i^ . Eq. 9.1.7, that it vanishes at the origin. This holds for real as well as imaginary argument. The behavior 
of the second Bessel function for small argument is given by 

Y.(^)«-^(|)"'r(z.), Rc(i.)>0, (34) 

cf. [i^, Eq. 9.1.9. Thus it diverges weaker than As the Bessel functions in (pk are multiplied by a factor 

this divergence is, however, cancelled. We thus conclude that the matter-dependent part of the wave 
function vanishes at the singularity. 

In the quantum mechanical case of regular potentials, one of the two independent solutions of the Schrodinger 
equation is not normalizable and therefore discarded. This is, however, not necessarily the case for vanishing 
angular momentum, since for the Coulomb potential, for example, both solutions are normalizable. The selection 
occurs in that case because only one of the solutions leads to an essential self-adjoint Hamiltonian [HO]- Such an 
argument can, however, not be invoked in the present case because the classical time parameter t is absent in 
quantum cosmology. There is thus no notion of unitarity here and thus no reason to advocate the self-adjointness 
of the Hamiltonian. We thus have non-uniqueness also in the "regular" case. 

• i - Va=0 < 0. 

This corresponds to the singular case in quantum mechanics. Here, vq is imaginary; we write i/q ~ ivq. Since 
the above formula ([55)1 holds for i/ ^ —1, —2, —3 . . ., we can still use it and find 

hm ^\3,,,X^l\m = lim ^ . e"^°'"i^j . (35) 

I^Ho I0HO r(l-|-u/o) 

For real as well as imaginary fc (corresponding to positive and negative energy, respectively) , this part of the 
matter-dependent wave function oscillates infinitely rapidly as it goes to zero as ]</>] — > 0. For the second part of 
the matter-dependent wave function we introduce a small parameter e > 0, since p4|) holds only for Re(i') > 0, 

hm hm v1^Y,+i,(]^]y£fc) = -- lim lim ^\ {^!^^ r(e + iz.o)e'"« '"('^) . (36) 

101^0 e^O TT |0|->Oe^O \ I j 
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This function oscillates very rapidly but goes to zero as \4>\ ^0 due to the factor y^j^f. We therefore conclude 
that the wave function vanishes at the origin in this parameter range as well. 

• i - Fa=0 = 0. 

In this case we have to consider the Bcssel functions Jo(z) and Yq{z) for small z. It is Jo(0) = 1 and so this 
part of Lpk vanishes when multiplied by -^/j^ffor 0-^0. But Yo(2:) w ^ In (z) for small z, [i^, 9.1.8, and thus 
diverges but ^/zYo{z) ~ ^y/z\n{z) — > as z — > 0. We thus find a vanishing wave function in this case as well. 

The singularities in these three cases are thus avoided when the matter-dependent part of the wave function vanishes. 
To complete this claim we have to insure that the gravitational part of the wave function docs not diverge at the 
respective singularities. Note that this result does not depend on any boundary conditions but is just a consequence 
of the Wheeler-DeWitt equation. 

Note also that the general solution ipk does not vanish at the singularity if v ^ i. This corresponds to = 0, thus 
a vanishing potential. This case is excluded by the set-up of our model. One could read this as: no potential implies 
no singularity and thus no singularity avoidance but a finite wave function at the origin = 0. Nevertheless, it is 
always possible to pick up a specific solution, that is, a particular (pk which vanishes. 

Concerning the gravitational part of the wave function, inserting the solution pi[l for the matter-dependent part 
of the wave function into the Wheeler-DeWitt equation ([24)) . we arrive at 

y {2Ck'pk + Ck(pk) + [^Ck + k^Ci^j y^fc = , (37) 

where a dot denotes derivative with respect to a. To obtain the gravitational part of the wave function, we assume 
that the terms Ck^k and Ck'fik can be neglected; this is just the meaning of the Born-Oppenheimer approximation 
discussed above. Basically, we assume that Ck varies much more rapidly with a than ipk and, moreover, neglect the 
back reaction of the matter part on the gravitational part. In summary, we are assuming that the change in the 
matter part does not influence the gravitational part; the matter part simply contributes its energy through k^. This 
leads to the equation 

(yC-fc + k^Ck^ ^k^O. (38) 

It is solved by 

C,(a) = 6ie'"?'" + 62e-'^" . (39) 

As the equation determining the gravitational part of the wave function is independent of £, we get the same solution, 
irrespective of whether we deal with a phantom or with a scalar field. Note that for < 0, k becomes imaginary and 
the dependence on a becomes exponential. In any case, Cfe(a = 0) < oo, so the wave function remains finite at the 
respective singularities and we can safely speak of singularity avoidance. 



V. BOUNDARY CONDITIONS AND SINGULARITY AVOIDANCE 



Nobody knows what the correct boundary condition for the quantum universe are. There have been several 
proposals, most of them using the boundary condition with the ambition to lead to singularity avoidance [H, [sil . 
|52|. |53||. DeWitt, in particular, speculates that the fact that a boundary condition is generally needed to make the 
quantum theory 'singularity-free' is an argument in favor of the theory: the theory itself does not leave any freedom 
of choice but provides the boundary condition itself [43| . As for the present state of understanding quantum gravity, 
however, singularity avoidance is demanded from the outside as a selection criterium. 

Irrespective of this discussion, it is in general necessary to impose a completely different condition if one instead 
wants to construct wave packets that follow classical trajectories with turning point in configuration space [i^ . 
Namely, one has to require that the wave packet decays in the classically forbidden region. This allows the interference 
of wave packets following the two branches of the classical solution behind the classical turning point. Of course, 
this is just the standard quantum mechanical treatment of classically forbidden regions. In general, out of solutions 
to the Wheeler-DeWitt equation which grow in the classically forbidden region, no wave packet can be constructed 
that follows the classical path. This is what happens generically for the no-boundary state (H^]. In order to make 
connection with the underlying classical theory, we shall therefore use herein as condition that the wave function 
decrease in the classically forbidden region. 
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One cannot assume that the condition of the wave function to decrease into the classicahy forbidden region always 
fixes it uniquely. Nonetheless it is a necessary condition that has to be imposed in order to get a theory with correct 
semiclassical limit. It can be used in selecting physically sensible solutions out of the general set of solutions to the 
Wheeler-DeWitt equation. We now want to explore the situation for our models presented above. 



A. Defrosting the big freeze 

1. Standard scalar field 

In the classical scalar field model, the region a < Omin = 1 is a forbidden region. We impose the boundary 
condition that the wave function decay there. In terms of the variable a, the boundary condition to be imposed is 
therefore ^ ^ as a — oo. But as a ^ —oo, — > i. The matter-dependent part of the wave function is then 
given by 

lim ipk{a,\(l)\) ^ J-[cisin{k\(j)\) - C2C0s{k\(f>\)] , (40) 

a — ^ — oo V TT 

see [i^, Eqs. 10.1.11 and 10.1.12. This vanishes for small |(?!)| if C2 =0 for both real and imaginary k. Thus imposing 
C2 = 0, we are left with the solution^^ 



(^fc(«,|0|) = ciV|0|J,(fc|0|) . (41) 

Due to the vanishing of tpk at the origin and its boundedness at infinity, the ipk are orthogonal, that is, J d4>ipkipn oc 
S{k ~ n) [5^. This is due to the fact that (j) can become arbitrary large. If the integration range of 4> was finite, 
we would need additional conditions to ensure orthogonality, cf. (49|, Eq. 11.4.5. Despite the orthogonality, k is not 
restricted to be of integer value. 

Inserting the solution (|1T|) for the matter-dependent part of the wave function into the Wheeler-DeWitt equation 
([24)) . we arrive at ((37|) . To obtain the gravitational part of the wave function, we proceed as before and find (|39p for 
the gravitational part of the wave function. 

For positive energy, k^ > 0, the gravitational part of the wave function is thus oscillating. No further complications 
arise and the full solution is given by 

*fe(a,</.) = civ^J,(fc|0|)[5ie'^" + 62e-'"^"] . (42) 

For negative energy, however, the gravitational part becomes exponential. To ensure that the boundary condition 
^' — > as a ^ —00 is satisfied for the entire wave function, we have to set 61 = 0. Thus, for imaginary k the 
gravitational part of the wave function decays exponentially for a — > —00 whereas the matter part remains finite, see 
Eq. (j4ip . The gravitational part alone ensures in this way that the wave function vanishes as a ^ —00. No additional 
condition arises for ipk- The full solution for imaginary k — > ifc is thus 

*fe(a,0) =62e"^"[ciJ,(ifc|0|) + C2Y,(ifc|(/)|)] , (43) 

and ipk and ipn are not orthogonal in this case. 

Notice that the functions ^'^(a, (j>) given in (|42)) and ([43]) fulfil as well the DeWitt criterium as they vanish for q = 
and 6 = 0. 



2. Phantom scalar field 

The classical model with the phantom-driven generalized Chaplygin gas has a classically forbidden region given 
by a > amax = 1. We therefore impose the boundary condition ^' — > as a ^ 00. In this limit, v becomes purely 



Recall that we use d = — and drop the tilde. 

In the quantum mechanics of a repulsive r~^-potential, r^/^J^(fcr) is the unique solution which vanishes at r = 0, cf. Equation (3.3) in 
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imaginary and large. We shall set here v := \v and v ^ oo (y now being real). We are thus looking for a combination 
of Bessel functions which vanish for large imaginary index. 

The wave function for positive energies, that is, /c^ > can be written as 

</7fe(a,|0|) = ^[ciJi,(ifc|,^|) + C2Yi,(ifc|0|)] . (44) 

As in this case the argument of the Bessel functions is purely imaginary, it is rather convenient to rewrite the wave 
function in terms of the linearly independent modified Bessel functions [49.] , 

¥.fc(a,|0|) = ^[ciIi.(fc|0|) + 52Ki,(fc|0|)] , (45) 

where ci, are arbitrary constants. It can be shown that the modified Bessel function Ki,^(fc|0|) goes asymptotically 
to zero as — > oo by using 



Ki, = V2(i 



X ) * exp \ — -^TT I X 



const. + 0(Sy^ - 2^^)"') 



1/ > x > 0, 



(46) 



cf. [5y], Eq. 7.13.2(19). Indeed, the modulus of the fimction Kiy(fc|(/)|) is oscillatory and its local extremum goes to 
zero as — !■ oo. Therefore, the implementation of the boundary condition gives 

(47) 



^>^{4>) = ^2^/\m^Ak\<\>\) ■ 

On the other hand, the wave function for negative energies, that is, < 0, can be written as 

^k{a, |(/)|) = ^ fciJi,(fc|0|) + C2Yi,(fc|,^|)l , 



(48) 



where k = ik and is positive. Here again it is more convenient to rewrite the general matter wave function in terms 
of Hankel functions in order to impose the boundary condition. We then rewrite the previous wave function as 



^,(a,|0|) - Vl-^l diii\'Jik\cl,\) + d2ll\:>{k\^\) 



t(2)/ 



(49) 



where di, d2 are arbitrary constants to be fixed by the boundary condition. It can be checked that Hj^^ (fc|0|) diverges 
for large v because 



TT 



exp 







"TT 


.71" 






exp 




1 — 


X 
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4J 





const. + £'((1/2+ a;2)-3) 



cf. [56|, equation 7.13.2(22). On the other hand, it can be shown that h|^' (fc|(/)|) vanishes at large values of 



v,x > 0, 
(50) 



This can be shown by combining the following properties of the Hankel functions: Il\'^\x) — cxp(— 7ri^)H'^/^(a;), 

ii-i^ix) = {H.ll}{x))*, where * denotes complex conjugation, and the asymptotic expansion of h|^^(x) for large order 
as shown above. Therefore, the wave function that fulfils the boundary condition for negative energies reads 

ifkiaM)=d2VW\il^\m)- (51) 

Before proceeding further, it is worthwhile to point out the following. The ordinary scalar field is defined on the 
entire real line. The phantom field, on the other hand, is restricted to the interval |(/)| G [0;(/)*]. Orthogonality does 
therefore not hold straight away but follows only if equation p2p is satisfied in the boundaries in which is defined. 
To obtain orthogonality, additional conditions would be necessary. As we cannot require that (pk vanishes at 0^,, 
as this would be in conflict with the requirement that the wave function follows the classical trajectory, we have to 
demand that its derivative with respect to (f> vanish there. This would require an analytic expression of the zeros of 
the first derivative of the modified Bessel function K^^ and the second Hankel function. To our knowledge, such an 
expression does not exist. Also, there is no physical motivation for this additional condition. 

The Cfc have to satisfy + k^Ck = 0. The solution is given by Eq. ([55)1 : (i) For positive energy, the imposition 

of the boundary condition leaves bi and 62 as arbitrary constants and therefore we need to impose the boundary 
condition on the matter-dependent part in this case, see Eq. (|47p . (ii) For negative energy, the imposition of the 
boundary condition picks up the exponentially decreasing function. The decay of the wave function for large a is 
thus already guaranteed through the purely gravitational part of the wave function. In principle, there is no need 
to impose the boundary condition on the matter-dependent part in this case, as long as the general solution for the 



r(2) 
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matter part is finite. However, it is not the case, so we have as well to impose the boundary condition on the matter 
part. This implies that only the second kind of Hankel function H^^'' is present on the matter part (see Eq. ([5T|l ). 
In summary, the physical solutions arc 

= (feie'^^^ + feae-'^") ^K,,(fc|0|), > 
*fc(a,0) = daexp l^-^fca^ V^H|^')(fc|0|), fc^ ^ 0. (52) 

These restrictions have to be taken into account when constructing wave packets (which will not be attempted here). 
Before concluding this subsection, we remark that the functions 'I'fc(a, (/)) given in ([5^ fulfil the DcWitt criterium as 
well, since they approach zero for a = and = 0. 



B. Smoothing the big-demarrage singularity 



As a third model with the same quantum structure, we consider the big-demarrage singularity. The equations for 
this singularity as it occurs for /3 ^ 1 are just the same as for the big-freeze generated by a phantom GCG presented 
above. The solution is also given by Eq. (|3ip . The difference lies in the classically forbidden region and thus in the 
boundary condition to be employed. Whereas in the previously studied phantom model we found a future singularity 
at a = flmax, this model has a past sudden singularity at a = Omin = 1. We therefore require the wave function to 
satisfy 'J ^ as a — > — oo. But as a — > — oo, i/ — > i. 

We split our analysis for positive and negative energies. For > 0, the boundary condition on the matter part 
of the wave function is Ci = 0, that is, we are left with the first Bessel function, while the gravitational part will be 
oscillatory. On the other hand, for fc^ < 0, the boundary condition on the gravitational part of wave function implies 
bi = 0, while no other boundary condition has to be imposed on the matter sector of the wave function. In summary, 
we obtain 



*,(a,0) = (6ie'"^" + 526-'"?^" ) ^\<l>\3.{tk\<j>\), 
'V6, 



fc^ > 



*fc(a,0) = exp — /ca ci J^(fc|0|) + C2Y^(fc|0|) 



< 0. 



(53) 



where k = ik. 

The solution is thus similar to the one found for the standard scalar GCG presented in the previous subsection. 
But note that is restricted to a finite interval here. 

Furthermore, as a consequence of the fact that the phantom scalar field has compact range, the matter dependent 
part of the wave function, ipk, is not orthogonal for different k. To obtain this, one has to require the condition (|32p . 

Here as well we notice that the functions ^'fc(Q:, 0) given in ([53]) fulfil the DeWitt criterium as they vanish at a = 
and 6 = 0. 



C. Other boundary conditions 

The Schrodinger equation with inverse square potential has been studied by various authors [3, [45| . They obtain 

as a solution the Hankel function of the first kind, H^,^^ for fc^ < and Va > \, ci. Equation (3.6) in Therefore, 
their solution is different from the ones we have obtained in the previous subsections. This is due to the boundary 
condition that is imposed in quantum mechanics: \1/ ^ as r ^ oo, r being the radial coordinate in the Schrodinger 
equation. In the quantum cosmological model this would correspond to the condition 5* — > as |(/)| — > oo. Whereas 
the vanishing of the wave function at infinity is a sensible requirement in quantum mechanics, it is not well- motivated 
in quantum cosmology. To demand a square-integrable wave function makes sense only in a Hilbert space with a 
probability interpretation and a (time) conserved probability. It is not obvious that we have such a structure in 
quantum cosmology: the Wheeler-DeWitt equation is timeless. 

Applying nevertheless this boundary condition to the quantum cosmological model with GCG mimicked by a 
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standard scalar field"'^^, we find (/3fc(a, (f) = ci •^/[^[Hi;^-' (ifc|0|) with = — fc^ and v purely imaginary as Va > |- This 
wave function corresponds to the choice C2 = ici. Requiring in addition orthogonality as in [4j], we would arrive at 
an energy spectrum given by 



En{a) 



2TTn 



■ exp 



n G Z 



(54) 



where kq is an arbitrary real constant. This can only be done for < 0. This is a similar non-uniqueness than in 
the quantum mechanical case [31 ■ 

Now, we can impose two different types of boundary condition on the total wave function '5: 



1. Cfc vanishes at the singularity; i.e. at a w 0, and (y3fc(a, 4>) is bounded at a w and (/> w 0, or 

2. Cfc vanishes well inside the forbidden region; i.e. a — > — oo, and (/?fc(a, 0) is bounded for a — 
We next impose both conditions separately. 



-oo and 



0. 



The equation for the gravitational part of general solution for the wave function obtained through a Born- 
Oppenheimcr approximation (sec Eq. p8|) ) can then be solved in the vicinity of the singularity, i.e. for a ~ 0. 
It satisfies 



Cfc - 7 (1 - 5a) Ck = 0, 



where the constants 7, 5 and c read 



6^2 



■ exp 



inn 



iTrnc 



(4c- 1)3 



Therefore, the gravitational part of the wave function for a w reads [4£ 



Cfe(a) = 61 Ai 



(^)^-(-57)*o 



6261 



(52 



(55) 



(56) 



(57) 



where Ai(z), Bi(z) denote Airy fmictions. Imposing that the wave function Ck vanishes at the singularity (a = 0) 
implies 



b2 = ~h — rf- (58) 

Bi[(3^)^_ 

It can be proven as well that (pk{o:, (p) = ci-\/[0[lli>^''(ifc|0|) vanishes for a w and 0^0. Therefore, we can conclude 
that the DeWitt criterium is compatible (in this case) with the boundary condition used in [3, ISl ■ 

The boundary condition s-Oasa— >— 00 cannot be imposed on the previous wave function (|57p . This is so 
because the wave function (|57p is valid only around a « 0. In order to impose the boundary condition 5" — > as 
a — *■ —00, we have to consider the wave function for a —00. In this limit, the energy spectrum (|54p reduces to^^ 



Ai 




Bi 





En{a) 



Consequently, 



Cfc(a) = hi exp 





+ 62 exp 


\/6ko \ 




a 









(59) 



(60) 



Notice that the boundary condition imposed in l44l l45f cannot be applied in our model with a phantom scalar field. The reason is that 
the scalar field is bounded in this case, so we cannot take the limit |(/)| ^ 00 in ipkip/-, <!')■ 

In taking this limit, notice that En{a) can become complex. However, En{a) given in II59I I is always positive. 
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Now, the condition Ck as a ^ — oo implies 62 = 0. On the other hand, ipk{o!, 4>) = ci yi0|Hi^^(ifc|0|) is bounded 
for a —00 and (/) 0. Therefore, we can conclude that the condition 4" — > as a — > —00 is fulfilled and compatible 
with the boundary condition used in [H, . 

In summary, as the wave function remains finite, the singularity would be avoided with the boundary conditions 
1. and 2. as well. But we emphasize again that the condition ^' — > as |0| — *■ 00 is perhaps not obligatory in the 
quantum cosmological case. It is presented here to compare our results to results presented elsewhere in the literature. 

VI. CONCLUSIONS AND OUTLOOK 

What are the main results of our paper? We have shown how the classical singularities of particular cosmological 
models describing dark energy features can be avoided in quantum cosmology. The framework is quantum geometro- 
dynamics, and restriction has been made to a class of models (containing a generalized Chaplygin gas) which could 
be of relevance for the description of dark energy. The wave function ^E* is defined on a two-dimensional configuration 
space consisting of the scale factor a and the scalar field representing the gas. After employing a Born-Oppenheimer 
type of approximation, we arrive at an effective stationary Schrodinger equation with singular (in the sense of quantum 
mechanics) potentials. The occurrence of such potentials in quantum mechanics signals an essential non-uniqueness 
of the solutions. The same happens here. Requiring that the wave function go to zero in the classically forbidden 
region still leaves a whole class of solutions which encode singularity avoidance in the sense of vanishing wave function 
(DeWitt-criterium) . This is different from our earlier papers [3 and [l^l in which we dealt with "regular" potentials 
and where appropriate normalizability assumptions singled out singular- free solutions (cf. also [Tojl. The case of the 
big-brake singularity in [23| contains a potential that behaves like the Coulomb potential in the vicinity of the classical 
singularity; as in quantum mechanics, no ambiguity remains after normalizability is being assumed. 

Our analysis has been performed in detail only for the case of the r~^-potential. A similar discussion can be made 
for the r~''^-potential leading to Mathieu functions; since the results are analogous to those of the r~^-potential, we 
do not present them here. We expect that a similar, although much more complicated, analysis can be performed for 
other singular potentials as well. It would be desirable to have a general classification of singularities with respect to 
singularity avoiding states at hand. 

We want to emphasize again one important conceptual aspect of this (and earlier) work. The singularities in these 
models do generally occur for a big universe, that is, for a value of the scale factor much bigger than the Planck scale. 
Singularity avoidance by quantum gravity then necessarily entails quantum effects for a big universe. Such quantum 
effects have been discussed earlier in the case of a classically rccoUapsing quantum universe [STj ; the quantum effects 
occur there near the (non-singular) region of the classical turning point as a consequence of the boundary conditions. 
We thus see that quantum cosmology may be of direct relevance for the future of our Universe. As discussed at length 
in [57} . the occurrence of such quantum effects possesses essential importance for understanding the arrow of time. It 
would be of great interest to study the question of entropy increase in the present models. 

There are many further open questions. In quantum mechanics, the remaining ambiguity of the solutions and the 
ensuing spectra points to new physics a short distances. The ambiguity must therefore be fixed either by reference to 
experimental results or by knowledge of the deeper theory. How is the situation in quantum cosmology? Experiments 
do not yet seem available, but can the ambiguity be fixed by theoretical means? One possible project would be to 
investigate whether prominent boundary conditions such as the no-boundary [sH or tunnelling (521 conditions are able 
to fix it. Perhaps the ambiguity can only be avoided by going to a more fundamental theory such as string theory. 
On the other hand, it may as well be imaginable that cosmological models leading to such singular matter potentials 
have to be excluded — an interesting selection critcrium for equations of state. 

An essential ingredient in our whole analysis is the demand that the wave function go to zero in classically for- 
bidden regions. How this relates in general to the hyperbolic structure of the Wheeler-DeWitt equation is not fully 
understood. Here we have avoided this problem by going to the Born-Oppcnheimer approximation. In general one 
would expect that certain quantization conditions on the physical parameters (masses, coupling constants, cosmo- 
logical constant) may appear [il], [sl, [s^, It would also be very interesting to establish a connection with a 
group-theoretical generalization of quantum geometrodynamics where it was shown that not only time, but even the 
ordinary three-dimensional space is absent at the most fundamental level, leading to singularity avoidance in a natural 
way [6l|. We hope that we can address some of these open issues in future publications. 

Finally, let us point out that a future sudden singularity has been recently investigated in (g^. The analysis was 
carried out in the loop quantum cosmology (LQC) framework^^. In this model, matter is modelled by a scalar field 



A discussion on the standard big bang singularity and LQC can be found in reference [3^ . 
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which rolls down (standard scalar field) or up (phantom scalar field) through an unbounded exponential potential. In 
the framework of LQC this scalar field behaves such that its effective energy density is finite; that is, there is a balance 
between its kinetic energy and its potential, while the effective pressure blows up in a finite future cosmic time. It 
was then concluded in [62l | that such a singularity cannot be avoided by means of the effective Friedmann equation in 
LQC. As it is mentioned in [6^ . the sudden singularity that appears for the phantom scalar field in the context of the 
modified Friedmann equation in LQC corresponds to a big-rip singularity in the standard relativistic case. For such 
a big-rip singularity, it was shown in [Tsj that it can be avoided in the sense that wave packets necessarily disperse 
when they approach the region of the classical big-rip singularity. It would be interesting to sec if this sort of big-rip 
singularity can be avoided as well by employing the DeWitt criterium that the wave function be zero at the classical 
singularity. 
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APPENDIX A: JUSTIFICATION OF THE GRAVITATIONAL WAVE FUNCTION APPROXIMATION 



In our models, we constructed the potential V{(j>) such that it corresponds to the polytropic equation of state for 
a GCG. The potential therefore depends on k and this dependence carries over to the matter-dependent part of the 
wave function. The Born-Oppenheimcr approximation can be understood as an expansion scheme with respect to k. 



The derivatives of tpk with respect to a are of non-zero order in k. This comes from tpk = I^i the vicinity 



. 2 

of the singularity, a « and — 
use to denote the value of at a : 



16 a-o Vt 

- 0. Therefore 



da 



Recall that Vi 



A 



^^^'^ and 



From Eq. ^ we see that Cfe - O (kO), Ck O [k-^) and - O (k-2) 
To obtain a consistent approximation scheme, such that in 



only the first term is kept and the others neglected, we have to assume that 



B I 3(1+3) 



Recall also that we 
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